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ABSTRACT 

A  general  superconvergence  result  of  conforming  finite  element  method  for  the  elliptic  problem  is  established  by 

using  L^-projection  method.  Regularity  assumption  for  the  elliptic  problem  with  regular  partitions  are  required.  Numerical 
experiment  are  given  to  verify  the  theoretical  results. 
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1.  INTRODUCTION 

The  finite  element  method  (FEM)  is  a  numerical  technique  for  solving  problems  which  are  described  by  partial 
differential  equation  or  can  be  formulated  as  functional  minimization.  To  summarize  in  general  terms  how  FEM  works  we 
list  main  steps  of  the  finite  element  solution  procedure  below.  The  first  step  is  to  divide  a  solution  region  in  to  finite 
elements.  The  finite  element  mesh  is  typically  generated  by  a  preprocess  or  program.  The  description  of  mesh  consists  of 
several  arrays  main  of  which  are  nodal  coordinates  and  element  connectivities.  The  second  step  selection  interpolation 
function.  The  interpolation  functions  are  used  to  interpolate  the  field  variables  over  the  element.  Often  polynomial  depends 
on  the  number  ofnodes  assigned  to  the  element.  The  third  step  is  to  find  the  element  properties.  The  matrix  equation  for  the 
finite  element  should  be  established  which  relates  the  nodal  values  of  the  unknown  function  to  other  parameters.  The  fourth 
step  to  find  the  global  equation  for  the  whole  solution  region  we  must  assemble  all  the  element  equations.  In  other  words 
we  must  combine  local  element  equations  for  all  elements  used  for  discretization.  Element  connectivities  are  used  for  the 
assembly  process.  Before  solution 

Boundary  condition  (which  are  not  accounted  in  element  equations)  should  be  imposed.  The  fifth  step  is  to  solve 
the  global  equation  system.  The  finite  element  global  equation  system  is  typically  sparse,  symmetric  and  positive  definite. 
The  sixth  step  is  to  compute  additional  results.  In  many  cases  we  need  to  calculate  additional  parameters,  for  example,  in 
mechanical  problems  strains  and  stresses  are  of  interest  in  addition  to  displacements,  which  are  obtained  after  solution  of 
the  global  equation  system. 

The  method  was  first  developed  in  1956  for  the  analysis  of  aircraft  structural  problems.  Thereafter,  within  a 
decade,  the  potentialities  of  the  method  for  the  solution  of  different  types  of  applied  science  and  engineering  problems 
were  recognized  over  the  years,  the  finite  element  technique  has  been  so  well  established  that  today  it  is  considered  to  be 
one  of  the  best  method  for  solving  a  wide  variety  of  practical  problems  efficiently.  In  fact,  the  method  has  become  one  of 
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the  active  research  areas  for  applied  mathematicians.  One  of  the  main  reasons  for  the  popularity  of  the  method  in  different 
fields  of  engineering  is  that  once  a  general  computer  program  is  written  it  can  be  used  for  the  solution  of  any  problem 
simply  by  changing  the  input  data.  Now  we  wall  talk  about  our  main  subject.  The  superconvergence  of  finite  element 
solutions  is  an  interesting  and  useful  phenomenon  in  scientific  computing  of  real  world  problems.  In  recent  years, 
supperconvergence  for  finite  element  solutions  has  been  an  active  research  area  in  numerical  analysis.  The  main  objective 
in  the  supperconvergence  study  is  to  improve  the  existing  approximation  accuracy  by  applying  certain  post  processing 

techniques  that  are  easy  to  implement.  The  superconvergence  is  obtained  by  applying  the  L^-projection  method  for  the 
finite  element  approximations  and  their  close  relatives.  In  this  article,  we  consider  the  poisson  problem  and  shall  develop  a 
general  superconvergence  results  for  it  is  conforming  finite  element  approximation. 

The  paper  is  organized  as  follows.  In  section  2,  we  present  review  for  the  elliptic  problem  and  the  method  to  find 
a  linear  system.  In  section  3,  we  develop  a  general  theory  of  superconvergene  result  by  following  the  idea  presented  in 
Wang  [3]. In  section  4,  several  numerical  experiments  are  present  to  support  the  theoretical  result. 


2.  PRELIMINARIES  AND  NOTATION 


Let  us  consider  the  poisson  problem  with  Robin  boundary  condition.  Let  Q.  be  bounded  and  open  subset  of  R 
with  Lipschitz  continuous  boundary (?fl. 


The  poisson  problem  seeks  an  known  function  'If  satisfying 


—Au  =  fin  n 


(2.1) 


— hu  =  gmdil 


(2.2) 


Where  T  and  ft  are  given  function.  Here  n  is  the  out  word  normal  vector  and"-  —        .  71  where        is  the 
gradient  of  u  and  _i  is  the  Laplacian  operator. 


Let 


space  of  square  integrable  function  define  on  £2  with  inner  product  and  norm 


l 


Further,  let  be  HillbertSobolev  space  of  order  7¥l,  that  is 


ff™(Q)  ={vEL2  (ft) 


DaV  G  L2  (Tt),  \a\  <  771  }With  inner  product 


i 


Iff  Ism  Q. 
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The  weak  formulation  of  (2.1)-(2.2)  is  to  seekUEH  ^(Tl)such  that 

a(u,  v)=  (f,  i?)n  +  {g,  vyaaVveH\tt)  (2.3) 

Where 

a  (u,  v)=  (Vu,  Vv)a  +  {u,  v]  m  =  fa  Vu.  Vvdx  +  J3a  u.  v  ds  (2.4) 

Let  Tfj  be  a  finite  element  partition  of  the  domain  £2  with  characteristic  mesh  size  h  and  lpr       be  the  set  of  all 
the  polynomials  defined  on  k,  with  degree  less  or  equal  to  r.  Define 

Vk  =  { veH1{ti)fv\k£P1{k}lVk€Th} 

The  finite  element  approximation  problem  for  (2.1)-(2.2)  is  to  find'lijj  £  Vfy  such  that 

a  (uhf  x?)=  (f,  v)a  +  is,  i?)9QyveVh  (2.5) 

Where 

a(uhl  v)  =  (Vuhl  Vv)a  +  (wfti        =  4  Vu*.  Fx?     +  /aQuh.  v  ds  (2.6) 

Note  that  Vfj  satisfies  the  following  approximation  property. 

^J^-^ll  +  /l^-^lll>^  W^Mh^M,  ^r+1(^)  (2-7) 
Define  a  norm  ||  -  ||  ]_  for  1/  —  /f  ^  (T?)  as  follows: 

1MB  =  a(>,  v)  =  Gfa,  Vx?)G  +  (y,  v)da  =  fa  Vv.  Vv  dx  -f-  /flQ  x?2rfs  (2.8) 

Definethe  basis  function  (p±  £  l1^  with  <J?j  =  Iwhere  i=j  and  (ff  j       ^  —  0  other  wise. 

Our  approximate  solution  lift  can  be  written  in  terms  of  it  is  expansion  coefficients  and  the  basis  function  as 
n 

i=1 

Expand  in  basis  lip,  insert  into  (2.5)  and  set  v=<Pi  ,  i  =  1,  ... ,  71 
n  n 

^  c,  (T^,  V^)n  +  ^  c(  C^ii  P/)«ki  =  if,  Vdo.  +  {9,  Vt)oa 

3=1  i=1 

This  is  a  liner  system  of  equations  AX=b  with  A=  (fly}  ,X=  fCj}  ,b=  {b^},foi  i  ,j=l,...,n 
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=  (f,<Pi)a  +  (S,<Pt)an 

Example  2.1:  Let  the  domain  O  =  [0,1]  X  [0,1]  .  Also  let  the  exact  solution  is  assumed  to  be 

w  =  smtjix)  cos(jiy) 

—Ait  =  2tt2  sinOrr)  cos  (Tiy)      in  ft 
du 

—  +  u  =  g     in  3D. 
on 

To  find  the  linear  system  Ax  =  b  see  Figure  1 


n3  < 


>  14 


n2 

Figure  1 

fly    =    CV^^)fl+  t<Pi,9jXi 

For  fej: 

^3  =  1  -  y,       then  V<p3  = 

For  &2  ■ 
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then  V<p4  =  [_^] 
then  SJ<p2  =  [  1] 


<p4  =  *  ~  y, 
<pz=l-x, 
<p2=y,      thenV<p2  =  [^] 
(Vqtlr  V<pO  =     Vft.  7<Pi  dfex  =  ^  .       .  = 
(V<p2f  V<p2)  =  /fcl  V<p2.V<p2  dkx  +  fk  ^  V<p2.  V<p2  dfc2 

i  rii  rii  .  1  roi  rQ' 


-iH-El-iGMSl- 


Then 


^1  = 


-1 

-1 

1 

0 

2 

2 

-1 

-1 

1 
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2 

-1 

-1 
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2 

-1 
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0 

7 

1 

5 

=  5lil-xfdx+£y2dy  =  \ 

l 


Then 
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A,= 


■  2  11 

366  ° 

12  1 

—  0  - 
63  6 

1  21 

6    °  36 

112 

0  

663 


Since  A  =  A±  +  A2 
Then 


r  s-i-i 

 o 

3  3  3 

—  15  -1 
 Cl   

3  3  3 

-1  S-l 

—  0  

3  3  3 

-1  -15 

0  

3    3  3  J 


A  = 


Now  to  find  the  load  vector 

„      \n  cos(nx)  cos  (ny)  1 
Vu=         .  \  \  .   '  A 
l—nsm{jix)  sin  Or}') J 

Then  /  =  2n2  sin(nx)cos  (ny) 


i  i 


b{  =  (J,  <pi)  =  j  /ipidfci  =  j  j  In2  sin  (tot)  cos(hj)  (y  —  x)dy  dx 


0  x 


_4 

7i 


b;  = 


f  f  -7T  -8+TT2 

=  (f,<Pi)  =  j  f92dk1+  J  f<p2dk2  =  —  +  - 


2n 
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b*  = 


71 

—K  -S  +  7T' 
2  27F 

2tt  2 
4 


Vit.  n2  +  u 


Vit.  n3  H-  u 


(— tc  sin  (rcy)  +  cos (Try)  sin  (Try) 
sin(?rr)  (irsm(rry)  +  cosijzy)) 
cos  (ttj)  (sin(Trx)  —  ncos(jzx)) 
cos  (Try)  (ttcos(tt?c)  +  sin(nx)) 


1  1 

=  j  -sin(nx) (1  —  x)  (te  +  J  —Ticos  (7ry}t£y 


Then  J&*  = 


1 
1 
1 


=  fcj       theji  x  =  A  1b 

Lemma  2.1 

There  exists  a  constant,  C,  independent  of  tl  such  th 

\a(vfw)\  <  CIMlillwlli,  Vw^H1^) 

Proof:  since  the  problem  (2.1)  -  (2.2)  satisfy  the  weak  form 
Where 
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1   Vw.  Vi?  dx 

+ 

1  w.vds 

a(w,v}  =  (Vw,  Vv)a  +  (w,  i?)3q  =  I    Vw.  Vt?  dx  +  I   iv.  i? 

By  using  Schwarz  inequality  we  obtain 

|a(w,t?)|  = 

^  l|Vw||i-(rolhMlj?"(H>  +  llwll^a^lMl^aQ) 

By  the  trace  inequality  there  is  a  constant,  C,  such  that 

IMIl2oq)  ^  c  IML 

We  have 

|a(w,v)|  <  IIVwIU^IIVvH^gd  +  C  HwlliClblli 

Since  1 1  VlV  ||Lz (Q>  <  llwllj 

l  l 
WhereUVwIL^  =  (ffl   \Vw\2dxf  ,\M\l=  (fn  \Vw\2  dx  +  f^ds)* 

Then  we  have 

k(w,iOI<  DwDJHIi  +  ^DwDilklli 

<L0-+C*)\M\i\\v\\i 
<  CllwII  Jblli 

Where     C  =  1  +  C2 

And  thus  the  proof  is  complete.  H 
Theorem  2.1 

Let  tt  and  lift,  be  the  solution  of  (2.1)- (2. 5)  respectively.  Then,  there  exist  a  constant,  Cf  independent  of 
h  such  that 

\\u-uh\\1^Chk\\v\\k+1_ 

Proof: 

Since 
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Hu-UftHi  <  inf  ||u-  f/lli 

and  by  using  (2.7)  then 

Hu-uJa  <  Chk+1\\u\\k+1forU€H^Hm 

3.  SUPERCONVERGENCE  BY  £2-PROJECTION 

-projection  technique  was  introduced  by  Wang  [3].  It  projects  the  approximate  solution  lift  to  another 
finite  element  dimensional  space  associated  with  a  coarse  mesh.  The  difference  in  size  of  the  two  meshes  can  be  used  to 

achieve  a  supperconvergence.  Now,  we  start  with  defining  a  coarse  mesh  7^  where  T  »  h  satisfying: 

T  =  ka  (3.1) 

With  CL  £  £0,1}.  For  now,  the  parameter  Ct  will  play  an  important  role  in  the  post  processing.  Define  a  finite 
element  space 

Herein  consists  of  piecewise  polynomials  of  degree  T  associated  with  the  partition  J^-. 

Let  Q?be  the  L^-projector  onto  finite  element  space V^-.    Ql^Can  be  considered  a  linear  operator  (projection) 
from  I?  (ft)  onto  the  finite  element  space  l^.We  will  assume  an  Hs,  S  ^  1,  regularity    for  the  solution  of 

(2.1),  (2.2): 

II  «  lls<  C  ||  /  |[s_2 

(3.2) 

Lemma  3.1 

Suppose  that  (3.2)  holds  with  1       S  5=  It  -+-  1  and  the  finite  element  space  1^  C  fis  ^         Then,  there 
exists  a  constant,  C7  independent  of  h  and  T,  such  that 

||  Q^U  -  QTUh  ||  <  Ch'-^^^WU-U^,  (3.3) 
Where  ^^(O.l^s  the  parameter  defined  in  (3.1). 

Proof:  The  definition  of  ||.  ||and  Q.z  give 

||  QTu  -  QTuh  ||=         SUp        ||  (Q^U-  QTUht  <p)|| 

From  definition  of  the  L^-projection  for  QT,  we  have 

(QTu -  Q.Tuhf  <p)  =  (u-  uhf  Q?<pl 
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Then 

II  Q^u-QTuk  ||=  snpvefi3(flj(l(pl=1||Cu-«fcj  QT<p)\\ ,  (3.4) 

Consider  the  following  problem:  find  W€fi  *  £flj  such  that 

a  (w,  v)  =  (QT  <p,  vyVvcH 1  (fi)  0.5) 

It's  follows  from  (3.2),wffs(fi)  f\  H  1  (fl)and 

||W||S<  lle^lL-z,  (3.6) 
For  some  constant C Using  the  inverse  inequality  on  the  right  hand  side  of  (3.6),  we  have. 

\\w\\s<CT^2-^\\<p\\o  (3.7) 
Subtracting  (2.5)  from  (2.3)  gives 

a{u-  Uklv)  =  0,  VV€Vh  (3.8) 
By  replacing  V  by  It      It in  (3.5)  and  using  (3.8),  we  obtain 

(u  -  uhl  QT<p)  =  a(u-  uhl  w) 

=  a{u  -  Uhl  W  —  1?)  (3.9) 
Where  V€Vh  .Using  Lemma  (2.1)  and  Theorem  (2.1)  we  obtain  from  (3.9) 

\(u-uh,QT<p)\  =  a\{u-uh,w-  v)\ 

<  C  inf  llu-UftllJIw-T?!!! 

^Ch^Wu-uMMU.  am) 

From  (3.6),  (3.7)  and  (3.1)  we  obtain 

\(u  -  uh,  QT<p)\  <  Ch?'1  ||  QT<p  ||s_2||  u  -  uh  ||x 
<  ck^H"11^0'2-^  \\u-uh  || -l ||  <p  || 

Then 

||  Gr^ft  II <  C^-1+minfe2-s}  ||  U-  IU.. 
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Theorem  3.1 

Assume  that  (3.2)  with  1  <  s  <  k  +  1  and  the  finite  element  space V^.  C         ^  (Tl}lf  the  exact  solution 
U  G  H k+1  (fi)  n  Hr+1  (Ti)  n  H  1(fl),  then  there  exists  a  constant,  Cf  such  that 

||u  -  Qxuh\\  4-  fc*|[Vr(u  -  QxUh)\\ 

<  Ck«(r+ti\\u\\r+1  +C^|t„-«ll||1  (3.10) 

Where  £7  =  S  ~  1  H"  Ct  111111  (0,2  —  S)  and  Uh  is  the  finite  element  approxima-tion  of  the  solutionis. 
Proof:  Since  Q^is  the  L^-projection,  then 

\\U-  QTU\\  *  CTr+1||u||rHhl  <  Cftff(r+l3||u||r+1.  (3.11) 
Since 

||it  -  QT«A||  <  ||«  -  QTu||  +  ||QT«  -  Qxuh\\ 

Combining  (3.3)  and  (3.11)  in  above  inequality  gives 

\\u-Qruh\\  <  Ch*(r+ti\\u\\r+1  +Ck'-x+aa**&2-*)\\u-viftif 

Which  complete  the  estimate  for  \\tL      QTUfa  ||  in  (3.10) 

Now  we  estimate  ||  VT  (VL  —  Q^lifr  || as  follows.  Since,  QT  is  L^-projection  and  using  (3.1)  we  have 
||VT(M  "  QxUh)\\  <  CTr||M||r+i  ^  Cft"-||tt||r+1.  (3.12) 
Multiplying  both  sides  of  (3.12)  by  flX ,  we  get 

h*\\VT(u-QTuh\\  <  £7i*<r+1)||ii||r+1  (3.13) 

The  inverse  inequality  gives: 

IIVr«?Tu -  Q7uh)\\  <  Ct^ICM* -  QTuh\\  (3  i4) 

Using  (3.1)  and  multiplying  (3.14)  by  h a the  above  inequality  implies 

ha\\Vx(Qvu-  QTuh)\\  <  C\\Qxu-QTuh\\ 

<  chs-1+<xmlR  <a2-*>  \\u  -  uJk,  (3.15) 

and 

h*\\VT{u  -  QTuJ||  <  fc"||VT(u- Qtm)||  +  A«|[V,  CQT«- QT«h)ll 
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Then 

II  («  -  QxuO  II        II  VT(u-  QTUjJ  || 

<  Cfoff(r+1]  ||  u  ||r+1+  C?ifc+s-1+amin^2-^||u||/.  +  1 

The  above  error  estimate  is  optimized  if  £t  is  selected  so  that 

a(r  +  1)  =  k+  s  —  1  +  a min(0,2  - s). 

Solving  CC  from  above  yields 
kA-s-1 

CL  =   ; — ;   (3.16) 

r  +  I— rain  (0,2— s) 

The  corresponding  error  estimate  is  given  by 

II  (u -  Qxu£  ||  +ha  ||  VT(u-  QTuh)  || 

<Ch^+1\\\u\\r+1+\\u\\k+1)  (3.17) 

Theorem  3.2 

Suppose  that  (3.2)  holds  with  1       $  +  1  and  the  finite  element  space  be  the 

finite  element  approximation  of  the  solution.  Then  the  post-processed  of  U  is  estimated  by  (3.17)  with 

k+s-1 


a  = 


r  +  1  -  min(0,2  —  s) 

4.  NUMERICAL  EXAMPLES  FOR  CONFORMING  FEM  BY  £f2-PROJECTION  METHOD 

In  this  section,  we  present  several  numerical  examples  to  verify  Theorem3.3.  The  triangulation  of  7\  is 
constructed  by:  (1)  dividing  the  domain  into  an  71    X  71   rectangular  mesh  and  (2)  connecting  the  diagonal  line  with  the 

positive  slop.  Denote  h,  =  —  as  the  mesh  size.  The  finite  element  spaces  are  defined  by 

■n3 

VT  =  {v€L2m:v\kePz(K),Vk6TT} 

Example  4.1:  Let  the  domain  fi  =  [0,1]  X  [0,1]  and  the  exact  solution  is  assumed  as 

u  =  (r3  —  x2')  (y3  —  y2') 

Table  1  shows  that  after  post-processing  method,  all  the  errors  are  reduced.  The  exact  solution  in  L^-norm  of 
\\u  —  Qx^h  Hi  has  the  similar  convergence 


Impact  Factor  (JCC):  4.2949 


Index  Copernicus  Value  (ICV):  3.0 


Super  convergence  of  Conforming  Finite  Element  Approximation  for  the  111 

I2 

Second  Order  Elliptic  Problem  With  Robin  Boundary  Condition  by  ir    -Projection  Methods 

Rate  as  ||u  —      ||.  There  is  no  improvement  for  the  U  in  L^-norm.  However,  the  error  in  the  M  ^-norm  has  a 

higher  convergence  rate,  which  is  shown  as  Q^tl^"^  for  ||VT  (u  —  QTZth)  ||,  see  Figure  2(a)  and  (b).  Figure  3(a)  and 
(b)  give  the  results  for  the  finite  element  approximation  of  (2.1)  and  (2.2)  before  and  after  post-processing. 


Table  1:  Error  on  Uniform  Triangular  Meshes  7},  And  TT 


h 

|  U  — Uk  |  i 

.  it  -  Uv  |  |K*  —  QTUh  |  -t  |  U  -  QTUh \_ 

V 

0.8936e-2 

0.6620e-3 

0.6455e-2 

0.6340e-3 

3"3 

0.2842e-2 

0.5986e-4 

0.1172e-2 

0.5707e-4 

4"3 

0.1212e-2 

0.1066e-4 

0.3510e-3 

0.1015e-4 

5"3 

0.6224e-3 

0.2795e-5 

0.1404e-3 

0.2661e-5 

6"3 

0.3604e-3 

0.9358e-6 

0.6698e-4 

0.8909e-6 

o(kklk= 

0.9752 

1.9915 

1.3873 

1.9933 

10" 
10": 
10- 

itf 

10"' 


-a-l|U-Uhl| 
-+-I|u-Quh|| 

»Q(h2) 


10" 


10" 


10  10 

(a) 


10° 


10"' 


-h-Iu-Qu^ 

v 

^O(h) 

10" 


10" 


10' 


10° 


(fe) 


(a)  (b) 
Figure  2:  (a)  Convergence  Rate  of  LJ-Norm  Error;  (b)  Convergence  Rate  of  ff  1-Norm  Error 


Figure  3:  Result  for  U  =  (x  3  -  X2}(y3  ~  y2)  In  Example  4.1.  (a)  Surface  Plot  of 
Approximation  Solution       (b)Surface  Plot  of  Approximation  Solution  Q^t^ 

Example  4.2:  Let  the  domain  ft  =  [0,1]  x  [0,1] 

and  the  exact  solution  is  assumed  as 

u  =  x5cos  (y) 

From  the  results  shown  in  Table  2,  it  is  clear  that  the  exact  solution  U  in  H  *-norm  has  the  superconvergence  but 
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there  is  no  improvement  in  L^-norm,  see  Figures  4  (a)  and  (b).  The  finite  element  solution  given  in  Figures  5  (a)  and 
(b).This  agree  well  with  the  theory. 

Table  2:  Error  on  Uniform  Triangular  Meshes      and  TT 
h  |  u  —  Uf..  |  1  |  it  —  Ui.,  |  |  u  —  QTuh  |  !  |  u  —  QTuh\\ 

3"30.3407e-10.5738e-3  0.8608e-20.4887e-3 
4"a0.1441e-10.1023e-3  0.2618e-20.8709e-4 
5~30.7384e-20.2684e-40.1065e-20.2288e-4 
6~a0.4274e-20.8991e-50.5135e-30.7648e-5 
0  (kk  ), k=  0.9982 1 .9986 1 .35561 .9993 


(a)    (b) 

Figure  4:  (a)  Convergence  Rate  of  Z/-Norm  Error;  (b)  Convergence  Rate  of  ff*-Norm  Error 


Figure  5:  Result  for  u  =  %5co-s  (y)  in  Example  4.1  (a)  Surface  Plot  of  Approximation 
Solution       (b)Surface  Plot  of  Approximation  Solution  QTUh 

REFERENCES 

1.  Nikishkov,  G.  P.,"Introduction  to  the  Finite  Element  Method",  Lecture  Notice,  UCLA,  (2001). 

2.  Singiresu,  S.  R.  The  Finite  Element  Method  in  Engineering,  Elsevier  Science  &  Technology  Books,  December 
2004. 

3.  Amiya,  K.  P.  Discontinuous  Galerkin  Methods  for  Second  Order  Elliptic  Boundary  value  Problems.  Department 
of  mathematics  Indian  Institute  of  Technology,  Bombay  Powai,  Mumbai -400076. 

4.  Per,  O.  P.  The  Finite  Element  Method-Lecture  Notes,  23  April  2013. 


Impact  Factor  (JCC):  4.2949 


Index  Copernicus  Value  (ICV):  3.0 


Super  convergence  of  Conforming  Finite  Element  Approximation  for  the  113 

I2 

Second  Order  Elliptic  Problem  With  Robin  Boundary  Condition  by  ir    -Projection  Methods 

5.  Hehu,  X.  and  Shaoqin.  G.  Superconvergence  of  the  Least-Squares  Mixed  Finite  Element  Approximations  for  the 

Second  Order  Elliptic  Problems.  2000  Mathematics  Subject  Classification.  65N30,  65N25,  31A30,  2  January 
2007 

6.  Wang,  J.  (2000)  A  super  convergence  analysis  for  finite  element  solution  by  the  least  squares  surface  fitting  on 
irregular  mesh  for  smooth  problems.  Journal  of  Mathematical  Study,  Vol.  33,  No.  3,  pp.  229-243. 

7.  Xiu,  Y."  Superconvergence  of  Nonconforming  Finite  Element  Method  for  the  Stokes  Equation",  Department  of 
Mathematics,  University  of  Arkansas  at  Little  Rock,  Little  Rock,  AR  72204-1099,  13  July  (2000). 


www.tjprc.org 


editor@tjprc.org 


